Percolation transition in networks with degree-degree correlation 
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We introduce an exponential random graph model for networks with a fixed degree distribution 
and with a tunable degree-degree correlation. We then investigate the nature of a percolation 
transition in the correlated network with the Poisson degree distribution. It is found that negative 
correlation is irrelevant in that the percolation transition in the disassortative network belongs to 
the same universality class of the uncorrelated network. Positive correlation turns out to be relevant. 
The percolation transition in the assortative network is characterized by the non-diverging mean 
size of finite clusters and power-law scalings of the density of the largest cluster and the cluster size 
distribution in the non-percolating phase as well as at the critical point. Our results suggest that 
the unusual type percolation transition in the growing network models reported recently may be 
inherited from the assortative degree-degree correlation. 

PACS numbers: 89.75.Hc, 64.60.-i, 05.70.Fh, 05.50.+q 
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I. INTRODUCTION 

Percolation in complex networks have been attracting 
a lot of interest in the statistical physics community [l|, 
A network may undergo a phase transition as one dis- 
cards nodes or links successively. When the fraction of 
remaining nodes or links is greater than a threshold, the 
network possesses a giant cluster which consists of a finite 
fraction of interconnected nodes. In the opposite case 
the giant cluster disappears and all nodes disintegrate 
into small clusters. It is called the percolation phase 
transition that separates the two phases. The percola- 
tion transition in complex networks, as well as in regular 
lattices pj , is interesting because of its relevance to ro- 
bustness of network systems against random failure and 
epidemic spreading [| S S S, 0, S S [M [H O • 

The random network of Erdos and Renyi (ER) is a 
prototypical model for complex networks (see Ref. 
for review). An ER network with N nodes is con- 
structed by linking each pair of nodes with the proba- 
bility p/[(N — l)/2], or by adding pN links between ran- 
domly selected pairs of nodes. The link density is given 
by p, and the degree distribution follows the Poisson dis- 
tribution Pdeg (k) = e~ ^ {k) k / k\ with the mean degree 
(k) = 2p. 

The ER network is uncorrelated in the sense that it 
lacks any structural correlation. This property allows 
one to study the percolation transition analytically. We 
summarize some known results: (i) The percolation order 
parameter P is defined as the probability that a node be- 
longs to a giant cluster. It exhibits a threshold behavior 
with the power-law scaling 
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(1) 



for p > p c = 1/2 with the exponent (3 = 1. (ii) Let 
n(s) be the number of clusters of size s per node. At the 
critical point it follows the power-law distribution 



with the exponent r = 5/2. For p =/= p c it does not follow 
the power law. (iii) The mean cluster size S is defined 
as the average size of finite clusters reached from nodes 
selected randomly. It also displays the power-law scaling 



S ~ \p - Pc 



(3) 



with the exponent 7=1. Note that the percolation 
transition belongs to the same universality class as the 
mean field transition [l|. 

The study has been extended to scale-free networks 
with the power-law degree distribution Pdeg.(k) ~ k~ x . 
Making use of the generating function method 0, [f| Q 
or the mapping to the q = limit of the q-state Potts 
model [1, [§] , researchers find that the percolation transi- 
tion in uncorrelated scale-free networks is characterized 
by the power-law scalings with the A-dependent expo- 
nents. 

Recently the percolation transition in a class of grow- 
ing networks draws interest [13l . [lj, I la . 1 1 €31 ] . The common 
feature of these networks is that the numbers of nodes 
and links are increasing in time with the density of links 
p fixed. Adding a node and making a link corresponds 
to nucleation of a cluster and merging of clusters, respec- 
tively. As one varies p, finite clusters condense into a gi - 
ant cluster giving rise to the percolation transition [17| . 
Interestingly the nature of the transition is different from 
that observed in the uncorrelated networks. The critical 
properties are summarized in the following: (i) The per- 
colation order parameter exhibits an essential singularity 



P ~ exp 



Vp-Pc 



(4) 



n(s) 



(2) 



with a constant a. (ii) The cluster size distribution fol- 
lows the power law n(s) ~ s~ r in the whole phase with 
p < p c - The exponent value is varying with p and takes 
t = 3 at the critical point with a logarithmic correction, 
(iii) The mean size of finite clusters S does not diverge at 
the critical point. Instead, it is finite and shows a discon- 
tinuous jump at p = p c . These features are reminiscent 
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of the Berezinskii-Kosterlitz-Thouless (BKT) transition 
in two dimensional equilibrium systems with continuous 
symmetry [3]. However, there is no similarity in the 
underlying mechanism for the transitions. 

Previous studies reveal that there exist at least two dis- 
tinct universality classes for the percolation transition in 
complex networks. One is characterized with the power- 
law singularity and the other with the essential singu- 
larity. It raises a question for the key ingredient that 
is responsible for the universality class. Similarly one 
may ask a question whether the essential singularity is 
observed in a non-growing network. 

There is an important observation that the grow- 
ing networks [HI, [3, EH, EH] have a positive degree- 
degree correlation. A positive degree-degree correlation, 
or phrased as an assortative mixing, refers to the ten- 
dency toward making links between nodes of similar de- 
grees [III . Consider a pair of connected nodes in a grow- 
ing network. As a network grows, the two nodes acquire 
more and more links generating a positive correlation. 
On the other hand, those networks displaying the power- 
law type percolation transition do not have any degree 
correlation. It gives a hint that the degree correlation 
may be an important factor determining the universal- 
ity class. In this work, we will investigate the effect of 
the degree correlation on the nature of the percolation 
transition. 

The degree correlation of a network can be quantified 
by the assortativity [l^| 

(fcfc'),-((fc + fc')/2) 2 
r ((k 2 + k' 2 )/2) l -((k + k')/2)y {0) 

where (•); denotes the average over all links and (k,k') 
denotes the degrees of two nodes at either end of links. 
Its sign indicates a positive (assortative) or negative (dis- 
assortative) degree correlation. It vanishes for uncorre- 
cted (neutral) networks. The degree correlation can also 
be monitored with the nearest neighbor deg ree AT nn (&) 
which is given by the average degree of neighbors of 
degree- fc nodes [2(| El- ^ is an increasing (decreas- 
ing) function of k for networks with a positive (neg- 
ative) correlation. Several characteristics of complex 
networks with the degree correlation have been stud- 
ied [H, [H, [24|, HH- However, the universality class of 
the percolation transition has not been understood yet. 

In this work, we will investigate the sole effect of the 
degree correlation on the nature of the percolation transi- 
tion. It necessitates a model for networks with a tunable 
degree correlation to a given fixed degree distribution. In 
order to avoid an interference with any other ingredient, 
the model is required to be random in other aspects than 
the degree distribution and the degree correlation. We 
propose such a model in Sec. |TTJ It belongs to a class of 
the exponential random graph model [2(|. In this class, 
a network model is defined as a Gibbsian ensemble of 
networks with an associated network Hamiltonian. The 
model and its structural properties will also be studied. 
In Sec. IIII1 we will investigate the percolation transition 



of the model as varying the degree correlation. Summary 
and discussions will be given in Sec. IIVI 

II. EXPONENTIAL RANDOM GRAPH MODEL 

The statistical ensemble approach is useful in modeling 
a network with a specific property [2(| [27], [H, [2i| . Sup- 
pose that one wants to construct a network model which 
is specified by an observable x. It is suggested that (26[ 
such a model can be defined as the Gibbsian ensemble 
over the set of networks Q — {G} with the probability 
distribution 

P(G) cx e- H ^ . (6) 

Here H(G), called the network Hamiltonian, is given by 

H{G) = Ox{G) . (7) 

The value of the observable x can be adjusted by the 
parameter 9 through the relation 

x=J2 P(G)x(G) . (8) 
Gea 

This is called the exponential random graph (ERG) 
model. 

Our purpose is to construct an ERG model for net- 
works with a fixed degree distribution Pdeg. (k) and with 
a tunable degree correlations. Then it may be natural to 
use the assortativity r in Eq. ([5]) for the network Hamil- 
tonian H. One can find a simpler form by using the 
requirement that the degree distribution Pdeg. (k) should 
be fixed. Note that ((fc + k')/2) l and ((fc 2 + k' 2 )/2)i are 
constants to a given degree distribution. Hence it suf- 
fices to consider the term (kk')i only in Eq. ([5]) for the 
Hamiltonian. 

Following is the formal definition of our model: Let Q' 
be the set of A-node networks that are specified by a de- 
gree distribution Pdeg.{k). A network G is conveniently 
described with the adjacency matrix A = (<%■) (i,j = 
1, • • • , A) whose matrix element takes the value ay = 
1 or if nodes i and j are connected or not. The model 
is defined as the Gibbsian ensemble over Q' with the net- 
work Hamiltonian given by 

7 N 

H{G) = --Y,<Hjkikj , (9) 

where ki = ■ ay denotes the degree of a node i and J 
is a control parameter. A positive (negative) correlation 
is favored by a positive (negative) value of J. The model 
may have any degree distribution. However, we only con- 
sider the simplest Poisson distribution as the ER network 
since we are interested in the effect of the degree corre- 
lation. 

The Gibbsian ensemble can be simulated by using a 
Monte Carlo method. We start with an ER network 



3 



FIG. 1: This picture shows a link rewiring process. Two links 
between node pairs (a, b) and (c, d) are chosen at random. 
They are then rewired to connect pairs (a, c) and (&, d). 




FIG. 3: Snapshots of networks with J = — 1 in (a), J = in 
(b), and J = 1 in (c). TV = 100 and po = 1. 
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FIG. 2: (a) Time evolution of the assortativity in networks 
with TV = 16000 (solid lines) and 32000 (symbols) nodes, 
(b) Stationary state values of the assortativity as a function 
of J in networks with TV = 32000 nodes. 



with TV nodes and L — poN links, and update net- 
work configurations via the so-called link rewiring pro- 
cess [3(| as illustrated in Fig. [TJ A link rewiring trial 
from a configuration G to G is accepted with the prob- 
ability min 

| l e -(ff(G')-ff(G))| The ^ the Monte Carlo 

dynamics leads to the Gibbsian ensemble in the station- 
ary state. It is noteworthy that the link rewiring process 
preserves the degree of each node. Therefore the ERG 
model combined with the Monte Carlo method allows us 
to study properties of complex networks with a given de- 
gree distribution but with a different degree correlation. 
The degree correlation can be adjusted with the param- 
eter J. 

Our model has a finite relaxation time. We tested re- 
laxation dynamics at N = 16000 and 32000 and po = 2. 
Figure [DJa) shows the time evolution of the assortativity 
r, averaged over N$ — 100 samples, at J = 1, 0, and — 1. 
One finds that the assortativity converges to stationary 
state values in finite Monte Carlo steps with a negligible 
finite size effect. 

The stationary state value of the assortativity at pa = 2 
is presented in Fig. EJb), which was measured with 
TV = 32000. We find that the assortativity vanishes at 
,/ = and is positive for J > and negative for J < 0. 
At J = 0, links are rewired randomly, which is supposed 
to lead to an uncorrelated network [30] - The assortativ- 
ity measure confirms the expectation. A positive (neg- 
ative) value of J leads to an assortative (disassortative) 




FIG. 4: The conditional probability distribution P(k'\k) for 
networks with TV = 10 4 and po — 2.0. J = — 1 in (a), J = 
in (b), and J = 1 in (c). 



network. 

Typical network configurations are shown in Fig. [3] 
They are obtained from the Monte Carlo simulations 
with J = — 1, 0, and 1, respectively, starting with the 
same initial ER network with TV = 100 and po = 1. 
Shown are only the largest cluster in each case. In the 
disassortative case (J = —1), most of large degree nodes 
with k > 2 are paired with small degree nodes with k = 1. 
On the contrary, in the assortative case (J = 1), large de- 
gree nodes and small degree nodes are segregated from 
each other. While large degree nodes form an interwoven 
core, small degree nodes form branches emanating from 
the core. The neutral network (J = 0) shows the features 
of the assortative and the disassortative networks simul- 
taneously. We note that the assortative network has the 
most inhomogeneous structure for the segregation. 

The degree correlation can also be seen from the proba- 
bility distribution p(k'\k) [l!|. It is the conditional prob- 
ability that a node at one end of a randomly chosen link 
has degree k' provided that a node at the other end has 
degree k. We find that it, being viewed as a function 
of k', is sharply peaked for all values of J. The peak 
position k' peak decreases, remains at a constant value, or 
increases when J < 0, J = 0, or J > 0, respectively. 
Numerical data showing these behaviors are presented in 
Fig. H 



III. PERCOLATION TRANSITION 

We proceed to study the percolation transition in our 
network model in the following manner: (i) An ER net- 
work is prepared with an initial link density po = 2. (ii) A 
correlated network out of the ER network is generated by 
applying the Monte Carlo dynamics to a given value of 
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FIG. 5: The density of the largest cluster P in (a) and the 
mean size of finite clusters 5* in (b) . These are obtained for the 
networks with N = 8 x 10 4 nodes averaged over Ns = 1000 
samples. 



J. (iii) Links are selected at random and removed suc- 
cessively. In the mean while, percolation properties such 
as the density of the largest cluster P, the average size 
of finite clusters S, and the cluster size distribution n(s) 
are measured as functions of remaining link density p. 
Those procedures are repeated Ns = O(10 3 ) times, and 
all measurements are averaged over those samples. 

We remark on the effect of the random link removal on 
degree correlation. Consider an arbitrary network with a 
link density po and an assortativity ro- Assume that the 
link density becomes p after the random link removal. A 
straightforward algebra shows that the assortativity r of 
the link-removed network is given by 



ro 



1 



l-p/po / 
p/po \ 



(fc 2 )/(fc)-l 



(10) 



where (k n ) is the n-th moment of the degree of the initial 
network [3l| . 

Our networks before the random link removal have the 
the Poisson degree distribution P d e a .(k) = e -^(k) k /kl 
with the mean degree (k) — 2po. So the moments are 
given by (fc 2 ) = (k) 2 + (k) and (fc 3 ) = (fc) 3 +3(fc) 2 + (ft). 
Inserting these into Eq. (ITU)) , one obtains that r = 
(p/po) r o- This relation guarantees that a network re- 
mains to be assortative (disassortative) during random 
link removals if it is assortative (disassortative) initially. 

We have studied numerically the percolation transition 
in the networks at several values of J. It seems that the 
nature of the percolation transition depends only on the 
sign of J. So, we will present the results for the cases with 
J = — 1, 0, and 1 as the representatives of disassortative, 
neutral, and assortative networks, respectively. 

In Fig. [SI we compare the density of the largest cluster 
P and the mean size of finite clusters S. For all values 
of J, the quantity P displays a threshold behavior indi- 
cating the percolation transition at a nonzero value of p. 
There are noticeable differences. The giant cluster shows 
up earliest in the assortative network (J = 1). It, how- 
ever, grows so slowly that it becomes smaller than those 



in the neutral (J — 0) and the disassortative (J = — 1) 
network at large p. These properties can be understood 
from the typical configurations given in Fig. [3] An assor- 
tative network consists of a highly interconnected core 
and branches emanating from it. The core is stable 
against the random link removal, whereas the branches 
can be easily disconnected from the core. Apart from the 
quantitative features, the data for P also suggest that the 
scaling behavior of P near the percolation threshold may 
be dependent on the assortativity. 

The behavior of S shows even more conspicuous differ- 
ence. There are sharp peaks near the percolation thresh- 
old for J = and J = —1. However, the assortative 
network with J = 1 does not exhibit such a peak. This is 
reminiscent of the percolation transition in the growing 
networks [3 Q EE El • 

Numerical data in Fig. 03a) and (b) indicate that the 
degree correlation may affect the nature of the perco- 
lation transition. We will investigate the nature of the 
percolation transition in each case using a finite size scal- 
ing (FSS) method. 

For finite values of N, the scaling law in Eq. (fT]) for P 
has the FSS form 

P(p,N) = N-V 7>((p-Pc)iV 1/P ) , (11) 

where v is the FSS exponent. The scaling function V{x) 
has the limiting behavior V(x ^> 1) ~ x@ and V(x — > 
0) = Ci with a constant c\. Similarly the scaling law in 
Eq. © for S has the FSS form 



(12) 



The scaling function S(x) has the limiting behavior 
Q(\x\ ^> 1) ~ x~ 7 and S(\x\ — > 0) = C2 with a constant 

C2- 

At J = 0, our model is equivalent to the ER random 
network. It is known that the percolation threshold is 
located at p c — 1/2. The critical exponents are given 
by those of the mean field theory, that is, f3 — Pmf = 
1 and 7 = jmf = 1 EL El- There is a little subtlety 
in the FSS exponent v. It was conjectured that v is 
given by the product of the mean field correlation length 
exponent vmf and the upper critical dimensionality d u 
provided that the criticality belongs to the mean field 
universality class [32|. The conjecture yields that v = 3, 
which is indeed the case for the ER random network 0] . 
In order to test the FSS ansatz, we have performed a 
scaling analysis. Figure [5] shows the scaling plots for P 
and S according to the FSS forms in Eqs. (TIT]) and (TT2"|) 
withp c = 1/2 and the mean field exponents /3 = 1,7= 1, 
and v = 3. All data taken from different network sizes 
N collapse onto single curves quite well indicating the 
validity of the FSS form and the critical exponents. 
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FIG. 6: Scaling plots for P and S at J = 0. 
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FIG. 7: Numerical results at J = — 1. (a) The mean size 
of finite clusters at various networks sizes N ranging from 
N = 10 4 (circles) to N = 64 x 10 5 (right-triangles), (b) 
The peak positions p pea k vs 1/iV ' 340 . The equation for the 
straight line is y = 0.553 + 0.628a;. (c) The peak height Smax 
vs. N. The straight line has the slope 0.340. 



A. Disassortative networks (J = —1) 

In this subsection, we will investigate the nature of the 
percolation transition in the disassortative network with 
J = — 1. In order to locate the percolation threshold p c , 
we focus on the FSS behavior of S plotted in Fig. EJa). 
It is evident that there are peaks which become sharp as 
N increases. If the percolation transition is character- 
ized by the power-law type singularity, the FSS form in 
Eq. (|12p implies that the peak position p max approaches 
the critical point p c as 



Pn 



Pi 



+ aN- 1 ^ 



(13) 



with a constant a and that the peak height S max grows 
as 



(14) 



We fitted the data S near the peak to a quadratic func- 
tion to interpolate the values of p ma x and S max at each 




(P-P C )N (P~P C ) N 

FIG. 8: Scaling plots for P in (a) and S in (b) at J = -1. 



value of N. Thus obtained values of p m ax and S max are 
fitted well to Eqs. (JT3J) and HU [see Figs. 7(b) and (c)], 
from which we find that 

p c = 0.553(5), 1/9 = 0.34(1), 7 /z7 = 0.34(1). (15) 

According to the FSS form in Eq. (fTTj) . one expects that 
the largest cluster P scales algebraically as 



at the critical point p = p c 
point, we obtained that 



N-P/» (16) 
. Fitting data near the critical 



P/v = 0.34(1) 



(17) 



Figure [8] shows that all data for P and S at different val- 
ues of TV collapse on to single curves, which proves the 
reliability of the numerical results for the critical expo- 
nents. 

The critical behaviors and the values of the critical ex- 
ponents are compatible with those of random networks. 
Therefore we conclude that the percolation transition in 
the disassortative network belongs to the same universal- 
ity class as that in the uncorrelated neutral network. 



B. Assortative networks (J = 1) 

In this subsection we will investigate the nature of the 
percolation transition in the assortative network. We 
have already noticed from Fig. [5] that the assortative 
network behaves differently from the neutral and the dis- 
assortative network. The difference is stressed again in 
Fig. [51 where we present numerical data for S obtained 
from networks at different sizes N = 10 4 ,- • • ,64 x 10 4 . 
Although there is a peak, it does not sharpen as N in- 
creases. At the same time, finite size effects arc non- 
negligible near the peak. This may be regarded as an 
indication that the assortative network does not undergo 
a percolation transition at all. It may be another possible 
scenario that there is a percolation transition associated 
with non-divergent S. 
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FIG. 9: The mean sizes of finite clusters in the assortative 
networks with J = 1 at different values of N. 
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FIG. 10: (a) P versus N at several values of p = 0.12,- ■ ■ ,0.24. 
(b) Effective exponent a versus l/N. 



We study FSS behaviors of the percolation order pa- 
rameter P. The FSS behaviors plotted in Fig. [TU] clearly 
shows that the network undergoes the percolation tran- 
sition at finite p c . As N increases, P(N) approaches a 
constant value for large p while it follows a power-law 
decay P(N) ~ N~ a for small p. We make use of an 
effective exponent a defined as 



a{N) 



\n[P{2N)/P{N)] 
hi2 



in order to locate the transition point. From the effec- 
tive exponent plot in Fig. [TUT b). we estimate that the 
transition point is at p c = 0.20(2). At the critical point, 
the density of the largest cluster follows the power-law 
scaling P ~ A~ Qc with 



a c = 0.6(1) 



(18) 



Note that the exponent a c is distinct from the corre- 
sponding value {(3/v) = 1/3 for the uncorrelated neutral 
network. Note also that S does not diverge at the perco- 
lation threshold. Based on these evidences, we conclude 
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FIG. 11: The cluster size distributions in the neutral net- 
works (J = 0) in (a) and the assortative networks (J = 1) in 
(b) of size TV = 10 6 . The straight line in (a) has the slope 
5/2, and the straight lines in (b) have the slope 3.9, 5.2, 7.1, 
respectively. 



that the percolation transition in the assortative network 
belongs to a distinct universality class. The assortativity 
is an essential ingredient for the universality class of the 
percolation transition. 

Our remaining task is to characterize the percolation 
transition in the assortative network. Figure [TUJ shows 
that P follows the power-law scaling P ~ N~ a not only 
at the critical point at p — p c but also in the non- 
percolating phase at p < p c . Furthermore the value of 
the exponent a is varying with p. It implies that the 
system is in a critical state for p < p c . The criticality is 
also observed in the power-law scaling of the cluster size 
distribution n(s) ~ s~ T with the continuously varying 
exponent r in the non-percolating phase. Figure [TTJ com- 
pares the cluster size distribution in the neutral network 
and the assortative network. In the neutral network, the 
cluster size distribution follows the power law only at the 
critical point with r = 5/2. On the other hand, it follows 
the power law both at and below the critical point in the 
assortative network. At the critical point, the exponent 
is given by r ~ 3.9. 

The power-law scaling behaviors of P and n(s) implies 
that the system is critical in the non-percolating phase. 
Hence the percolation transition cannot be described by 
power-law type scaling laws. Instead, the assortative net- 
work model shares many features in common with the 
growing network models [HI, [TU, [TU, [l6|] in regard to the 
critical behaviors. The non-divergence of S at the crit- 
ical point and the power-law scaling of P ~ N~ a and 
n(s) ~ s~ r in the non-percolating phase are such com- 
mon features. At the critical point our numerical esti- 
mates are a ~ 0.6 and r ~ 3.9, while the corresponding 
values are a = 1/2 and r = 3 in the growing network 
model [HI] . We attribute these discrepancies to the loga- 
rithmic corrections at the critical point [l5| . Our model 
is a generic one for networks with assortative degree cor- 
relation. Therefore our numerical results suggest that 
assortative degree correlation is responsible for the un- 
usual scaling behaviors observed in the growing network 
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models. 



IV. SUMMARY AND DISCUSSION 

In summary, we have investigated numerically the na- 
ture of the percolation transition in networks with de- 
gree correlation. As a model for the correlated net- 
works, we have introduced the exponential random graph 
model with the Hamiltonian given in Eq. ([9]) under the 
restriction that the degree distribution is fixed. Using 
the model combined with the Monte Carlo method ex- 
plained in Sec. [TTJ, one can generate correlated networks 
to a given degree distribution (taken as the Poisson dis- 
tribution in this work). Numerical results show that the 
negative degree correlation is irrelevant in that the dis- 
assortative network exhibits the same type percolation 
transition as the neutral network. On the other hand 
the positive correlation turns out to be relevant. The 
percolation transition in the assortative network is char- 
acterized by the non-diverging S at p = p c and power-law 
scaling of P ~ N~ a and n(s) ~ s~ T with the continu- 
ously varying exponents a and r in the non-percolating 
phase. 

The scaling behaviors of the assortative network are 
compatible with those of the growing network mod- 



els [ill, 03, [ID, [l6| . It strongly suggests that the unusual 
percolation transition in the growing network models is 
inherited from the assortative degree correlation. This 
conclusion is highly plausible but not decisive yet. It is 
worthwhile to mention a discrepancy in the property of 
the mean size of finite clusters S. The growing network 
models show a discontinuous jump in S at p = p c . How- 
ever we do not find an indication of such a discontinuity 
in the assortative network. It remains as an unsolved 
question whether the discontinuous jump in S is an uni- 
versal property or not. Numerical studies were limited 
to networks up to size N = 64 x 10 4 because the Monte 
Carlo dynamic generating correlated networks is slow. 
Numerical data up to that size fail to justify exclusively 
the essential singularity in P as in Eq. Q. In this re- 
spect, it is desirable to find an efficient algorithm with 
which one can generate the correlated network of much 
larger sizes. At the same time, it will help us understand 
better the property of correlated networks if one can find 
an analytically tractable model. These are left for future 
studies. 
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